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We have studied graphene/gold interface by means of density functional theory (DFT) and scan- 
ning tunneling spectroscopy (STS). Weak interaction between graphene and the underlying gold 
surface leaves unperturbed Dirac cones in the band-structure, but they can be shifted with respect 
to the Fermi level of the whole system, which results in effective doping of graphene. DFT cal- 
culations revealed that the interface is extremely sensitive to the adsorption distance and to the 
structure of metal's surface, in particular strong variation in doping can be attributed to the specific 
rearrangments of substrate's atoms, such as change of the crystallographic orientation, relaxation 
or other modifications of the surface. On the other hand, STS experiments have shown the presence 
of energetic heterogeneity in terms of the changes in the local density of states (LDOS) measured at 
different places on the sample. Randomly repeated regions of zero-doping and p-type doping have 
been identified from parabolic shape characteristics and from well defined Dirac points, respectively. 
The doping domains of graphene on gold seem to be related to the presence of various types of the 
surface structure accross the sample. DFT simulations for graphene interacting with Au have shown 
large differences in doping induced by considered structures of substrate, in agreement with experi- 
mental findings. All these results demonstrate the possibility of engineering the electronic properties 
of graphene, especially tuning the doping across one flake which can be useful for applications of 
graphene in electronic devices. 

PACS numbers: 73.22.Pr, 71. 15. Mb, 73.20.Hb 



I. INTRODUCTION 

Graphene 1 , a single layer of carbon atoms arranged 
as those in one sheet of graphite, has shown fascinat- 
ing physical properties 2 - that are currently being inves- 
tigated in basic research and which make it a promising 
material to play a key role in future microelectronics, op- 
toelectronics and sensors^. Realization of most of the 
industrial applications requires reliable methods for pro- 
duction of large- area graphene, which has been practi- 
cally achieved via, for example, chemical vapor deposi- 
tion (CVD) growth on transistion metals^ — as well as by 
ethylene irradiation in the case of graphene on copper and 
gold substrates^. However, the interaction with polycrys- 
talline metallic foils can strongly affect the graphene's 
electronic characteristics and impair the performance of 
potential devices. Similarly, the presence of metallic con- 
tacts which are essential in electronics also influences its 
fundamental properties. Due to all these facts, mecha- 
nisms of coupling between graphene and metals is now 
one of the leading fields in the research of graphene. 

DFT calculations and experimental studies by angle- 
resolved photoemission spectroscopy (ARPES) have 
shown that graphene can bind strongly or weakly on var- 
ious metallic substrates. One can distinguish two main 
types of interaction between the graphene sheet and the 
surface of a metal, i.e. chemisorption that completely 
eliminates the characteristic Dirac cones from the band- 
structure, as for example in the case of Co, Ni and Pd 
substrates^ — , and physisorption which leaves the linear 
dispersion undisturbed, as observed for Al, Ag, Cu, Au 



and Pt (111) surfaces. In the case of physisorption, the 
weak interaction causes the Fermi level shift downwards 
(upwards) with respect to the Dirac point. It results in 
doping by the holes (electrons) transferred from a metal 
to graphene which becomes p-type (n-type) doped. Ac- 
cording to most of the theoretical studies^ 2 - — , graphene 
on Al, Ag and Cu is n-type doped, while interaction with 
Au and Pt(lll) surfaces seems to cause p-type doping. 
Although the particular DFT calculations can lead to 
different values of the doping level for the same metal, 
it is usually assumed that one type of substrate implies 
definite type and level of doping and any contradictions 
are related to the limitations of ab initio methods whose 
results depend on the choice of the exchange-correlation 
interaction and on the calculation strategy^ 7 -. 

In our recent work, we have analyzed the correlation 
between the adsorption geometry and electronic prop- 
erties of the graphene/Au(lll) interface and we have 
shown the influence of the geometry parameters char- 
acterizing the modeled system 17 . We focused mainly on 
the aspects related to the strategy of the calculations, i.e. 
it has been found that the factors such as the in-plane 
cell parameter and interlayer spacing in gold can impact 
the Fermi-level shift and even a small change of geome- 
try may lead to a transition between p-type and n-type 
doping. Furthermore, the adsorption distance which has 
the strongest influence on doping level (see Fig. 3 in Ref. 
Il7l ) can be, in practice, set arbitrarily, since the equilib- 
rium position of the graphene sheet with respect to Au 
depends directly on the choice of one from many possible 
exchange-correlation functionals. 
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However, it should be stressed that tuning the geom- 
etry details may be realized experimentally and should 
then lead to predicted abrupt modifications of electronic 
properties. For example the presence of steps, buckling or 
corrugation of the surface could induce local variations 
in graphene-substrate distance. In this case graphene 
can be suspended above certain regions and coupled in 
a different way to the metal. Next, the change of the 
substrate's in-plane lattice parameter can be achieved 
by depositing graphene on another crystallographic face 
which is characterized not only by different lateral ar- 
rangements of metal atoms but also by different size of 
the unit cell of the interface. Moreover, the strong relax- 
ation of the surface can be the consequence of interaction 
with impurities present on the surface or may depend on 
details of the sample preparation process. 

All these surprising theoretical results motivated us to 
experimental studies of graphene/gold system in order to 
verify the extreme sensitivity of the doping level on the 
interface structure. Since the variations of doping are 
expected to be local, the most suitable tool for these in- 
vestigations are scanning probe methods. In particular, 
the STS mode provides information about first derivative 
of the tunneling current with respect to the bias voltage 
which is the measure of the local density of states at ev- 
ery point above the surface. The electronic structure is 
recorded locally, thus, in contrast to ARPES technique, 
one can detect nonhomogeneous electronic properties at 
different places on the sample and identify repeated do- 
mains of doping. 

Indeed, the preliminary STS measurements performed 
on graphene/gold sample have shown very strong ener- 
getic heterogeneity in terms of changes in LDOS recorded 
at different places on the surface 2 ^. Now, we have ob- 
served separate regions with completely different elec- 
tronic structures as well as slight changes in LDOS within 
single domains. The heterogeneity of electronic proper- 
ties seems to be domains of doping, but: (i) differen- 
tial tunneling spectra are not sufficient to give straight- 
forward information about doping level unless they are 
completed by another experimental or theoretical tech- 
nique which enables identification of the Dirac points in 
the profiles (see Refill! for wider discussion), (ii) there is 
no evidence which particular features in the structure of 
surface lead to the given changes in the doping level. It 
seems that the domains of doping can be created sponta- 
neously on metallic surfaces, but more knowledge about 
mechanisms of interaction is needed to control the doping 
or to create multi-doped samples for particular purposes. 

In the present paper, we attempt to make a link 
between theoretical models studied by means of DFT 
method and STS results obtained for graphene/gold sys- 
tem. First, we have thoroughly analyzed STS profiles of 
graphene/ Au system and determined the character and 
distribution of domains with different electronic proper- 
ties. Secondly, we have recorded dl/dV profiles in the re- 
gions on the sample which were not covered by graphene. 
This provides information of the structure of the pristine 



Au surface that needs to be included in the theoretical 
models. It should be stressed that the role of DFT simu- 
lations is twofold: (i) to identify unambiguously the Dirac 
points in experimental STS profiles which is achieved 
by providing simultaneously the information about the 
band-structure and LDOS above the surface according 
to the scheme proposed in Ref. 18, (ii) to show the 
electronic properties of graphene deposited on structures 
constructed in such way that the properties detected by 
STS for gold substrate were accurately described. Fi- 
nally, both simulated and experimental data are com- 
pared and assessed to be in a satisfactory agreement. 

The paper is organized as follows. In Sec. II the calcu- 
lation method and experimental details are summarized. 
In Sec. Ill we report on the results of STS measurements 
of the graphene/gold sample and on the dl/dV profiles 
recorded for pristine surface of gold. The theoretical in- 
terpretation of the observed heterogeneity is proposed, 
analyzed and compared with experimental data in Sec 
IV. Some final remarks and perspectives are discussed in 
Sec. V. 



II. METHODS 
A. DFT calculations 

First principles calculations based on the DFT have 
been performed using VASP softwar e - 21 ! 22 within the gen- 
eralized gradient approximation in the parametrization of 
Perdew, Burke and Ernzerhof (PBE) 23 empirically cor- 
rected to include long-range dispersion according to a 
scheme proposed by Grimm o 24 i 25 . The dispersion coef- 
ficient Cq for gold has not been included in the original 
paper of Grimme— , thus we have used a value of 40.62 
Jnm 6 /mol and of 1.772 A for the vdW radius of Au (iio), 
as proposed in Ref. [26|. The pair interactions up to a ra- 
dius of 12 A have been included in the simulations and 
the global scaling factor sq has been set to 0.75 since the 
PBE exchange-correlation was employed. Moreover, the 
projector augmented wave (PAW) pseudopotentia b 27 i 28 
were used for electron-ion interactions and the electronic 
wave functions have been expanded in a plane-wave ba- 
sis set of 400 eV. The electronic self-consistency criterion 
has been set to 10 _T eV. 

The (111) surface of gold is modeled by a periodic slab 
geometry. Supercells consisting of six/seven layers of Au 
atoms and a single sheet of graphene are separated by the 
vacuum thickness of at least 15 A placed in the direction 
normal to the graphene's plane to avoid interactions with 
spurious images of the slab. We employed the typical in- 
plane adsorption geometry as shown in FigHJ where 2x2 
graphene's and substrate's R30°(\/3 x y/3) unit cells are 
directly matched. The size of the supercell is determined 
by the in-plane lattice parameter ah ex denoted in FigJTJ 
We have performed calculations only for one value of «h e x 
adapted to the one of graphene which was optimized us- 
ing DFT-D2 method and estimated to be 4.932 Ail. The 



3 




Au 



(b) 

Au 



MG 




FIG. 1: (Color online) Top view (a) and side view (b) of 
adsorption geometry of graphene on a Au(lll). Carbon 
atoms are denoted as red (darker) balls, gold atoms as yellow 
(lighter) ones. Parallelogram defines the unit cell. 

structural relaxations have been performed for the un- 
covered substrate first, next the graphene sheet has been 
added above, and all atoms in the top two layers of metal 
as well as all carbon atoms have been allowed to move. 
Total energies were converged to within 10 ~ 6 with respect 
to the ionic steps. The determined equilibrium distance 
is equal to <i eq =3.23A. Other models of graphene/ Au in- 
terface are considered and described in details in Sec. 
IV. 

The simulations of the LDOS of graphene/ Au system 
have been done according to the scheme proposed in 
Ref. 18. The typical STS profiles in Tersoff-Hamann 
approach are obtained by an evaluation of the constant- 
height charge images for several values of the bias voltage 
followed by their numerical differentiation. It allows us 
to calculate the spectrum in the whole space above the 
crystal. The choice of the particular point above the sur- 
face corresponds to the STS data taken under open loop 
conditions at the fixed position. The spectroscopy sim- 
ulations for graphene on metallic substrate are usually 
demanding: at least one k-point every few meV for a 
specific band is needed^ which means that a very dense 
k-points mesh must be used in the final calculations of 
density of states. We use the tetrahedron scheme^ and 
the T-centered 24x24x1 k-point mesh during the opti- 
mization and self-consistent runs for accurate Brillouin 
zone integrations. 



B. Experimental details 

The STM/STS studies of monolayer graphene on poly- 
crystalline gold have been done on the sample prepared 
by Graphene Industries, UK. The sample consists of 
graphene flakes transferred on 8 nm Au with 0.5 nm 
Cr adhesion layer sputtered onto 100 nm Si02 grown 
in Si(001). The 8 nm gold layer is thick enough to 
be stable during graphene flakes deposition and suffi- 
ciently thin to be transparent for the visible light. Thus, 
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FIG. 2: (Color online) (a) Optical image of graphene flakes 
deposited on Au/Cr/Si02/Si substrate, (b) 450 nm x 450 nm 
STM topography showing the details of MG and Au border 
line, (c) dl/dV (E, line) map recorded on MG/Au interface 
along arrow shown in figure (b). Colors: blue, green, red - 
low, intermediate, and high value of the LDOS, respectively, 
(d) The example of dl/dV profile recorded on MG/Au sample. 



it is possible to extend the three-interface Fresnel-law- 
based model effective for optical graphene identification 
on widely used air/MG/Si02/Si samples^ to the five- 
interface system air/MG/Au/Cr/Si0 2 /Si^. This en- 
ables the optical identification of graphene flakes after 
deposition on polycrystalline gold. All the STM/STS 
experiments have been carried out at room temperature 
in UHV condition using VT-STM/AFM microscope inte- 
grated with the XPS/UPS/AES/LEED/MULTIPROBE 
P system (Omicron GmbH). The tips have been pre- 
pared by mechanical cutting from the 90%Pt-10%Ir al- 
loy wires (Goodfellow) and electrochemical etching of the 
W wires (Goodfellow). In the STS mode the I/V curves 
were recorded simultaneously with a constant current im- 
age by the use of an interrupted- feedback- loop technique. 
The obtained data have been used to calculate the first 
derivative of the tunneling current with respect to the 
voltage (dl/dV) and to build spatial conductance maps 
i.e. dI/dV(x,y,E). 



III. EXPERIMENT: RESULTS AND 
DISCUSSION 

A. STM/STS of monolayer graphene deposited on 

gold 

The monolayer graphene (MG) deposited on gold sub- 
strate has been initially identified by the visual inspec- 
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tion with an optical microscope. As it was mentioned 
above, the gold layer does not obscure the optical inter- 
ference effect that makes graphene visible (see Fig|2](a)). 
Next, the STM topography has been recorded at sample 
bias U=+0.8 V and with the tunneling current set point 
equal to 0.2 nA. The 450 nm x 450 nm topography is 
presented in Fig|2jb) which clearly shows the border be- 
tween MG/Au and pure Au. It seems that the graphene 
layer does not map exactly the structure of Au substrate. 
The estimated value of r.m.s. calculated over 200 nm x 
200 nm area on MG equals 0.5 nm while a similar mea- 
surement on gold gives value equal to approximately 0.86 
nm. 

In Fig|2](c) we present dl/dV (E, line) map calculated 
as a function of the energy (y axis) and of the position 
along selected line indicated by the arrow in Fig|2jb) (x 
axis). The most visible feature in the plot is the bor- 
der between MG/Au and pure Au regions. It is denoted 
in Fig(2fc) by vertical white line dividing the plot into 
two parts: the left-hand one shows MG deposited on Au, 
while the right-hand one is associated with uncovered 
Au substrate. It is clear that the major change in LDOS 
starts to appear at MG/Au - Au border. Any strong 
changes cannot be observed at the occupied part of the 
spectra, while the dramatic change can be easily noticed 
in the unoccupied part of the profile, where a considerable 
decrease in LDOS occurs at about 0.5 eV and is accom- 
panied with the presence of local maximum located at 
energy equal to 0.37 eV above the Fermi level. The rep- 
resentative dl/dV curve recorded in this region shown in 
FigEfd) illustrates the asymmetry of the profile with re- 
spect to the Ef. The minimum at 0.5 eV seems to mark 
the position of the Dirac point (Ed) of graphene, but 
only the comparison with DFT results allows for its pre- 
cise identification^. Moreover, the energetic position of 
Ed varies across the sample. Taking into account 65536 
individual dl/dV curves we estimated that the position 
of Ed is located in the range of 0.25 - 0.55 eV above the 
Fermi level. 



B. Heterogeneity of graphene's electronic structure 

The STM topography shown in Fig|2] (b) suggests 
strong stuctural heterogeneity of the substrate, in par- 
ticular the presence of numerous tiny crystallites of Au, 
which can be easily observed in both regions, covered and 
uncovered. Since the grains of Au may differ from each 
other, it can lead to the diversity in the local electronic 
properties. We have analyzed the considered MG/Au 
system using the conductance maps measured across the 
sample. In Fig|3]we present 100 nmxlOO nm dl/dV (x, 
y, E) maps at selected energies around the Fermi level. 
It is clear that the spatial distribution of the LDOS am- 
plitudes recorded at the same energy is not uniform and 
varies from region to region. It leads to the conclusion 
that MG is not homogenous in terms of local electronic 
structure, especially in the case of energies very close to 
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FIG. 3: (Color online) 100 nm x 100 nm dl/dV (E, x, y) 
maps recorded on MG/Au. Colors: blue, green, red - low, 
intermediate, and high value of the LDOS, respectively. 
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FIG. 4: (Color online) (a, b) dl/dV profiles recorded in differ- 
ent places on MG/Au sample. The dashed green lines corre- 
spond to the points in regions denoted as blue in Fig[3j while 
the spectra taken in the points within orange/yellow regions 
are marked as solid red ones. 



the Fermi level. At energies higher than Ep ±0.250 eV 
spatial distribution of the LDOS seems to be much more 
uniform and it is difficult to distinguish between differ- 
ent domains. The dl/dV profiles measured in the regions 
with high LDOS very close to the Fermi level (yellow 
and light blue patches in Fig. [3j) show well defined Ed 
located in the range of 0.25 - 0.55 eV above the Fermi 
level (see solid red lines in Figs|4ja),(b)). We have esti- 
mated that the size of domains with well defined max/Ex> 
feature in conductance maps typically spread over 10-25 
nm. The dl/dV curves recorded in the dark blue re- 
gions (low LDOS close to Ef) strongly differs from the 
profile presented in Fig. [2] (d), i.e. parabolic shape of 
dl/dV and the lack of well defined Ed point is observed 
(green dashed lines in Figs|4](a), (b))). It should also be 
stressed that, surprisingly, the existence of considerable 
LDOS heterogeneity close to the Fermi level allows us to 
easily identify regions on MG/Au/Si02/Si surface with 
well defined Ed, even though the Dirac point is located 
at about 0.25-0.55 eV above Ef- It can be explained by 
the fact that the presence of Ed is associated with high 
value of LDOS at the Fermi- level (see red solid lines in 
Figs. Ha),(b)). 

It can be easily noticed in Fig. |3]that apart from large 
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FIG. 5: (Color online) (a, b, c, d) dl/dV (E, line) maps 
recorded on different MG/Au regions. Colors: blue, green, 
red - low, intermediate, and high value of the LDOS, respec- 
tively. 



domains one can distinguish yellow/orange and light blue 
patches. Figures 2] (a), (b) show two different dl/dV pro- 
files (red solid lines) taken in two selected points that 
correspond to the additional fine structure visible in the 
maps presented in Fig|3j These two spectra presented in 
Fig. 2] (red solid lines) differ from each other even in the 
vicinity of the Fermi level, which is especially pronounced 
while compared with dl/dV curves (green dashed lines in 
Fig. [4j) taken in regions denoted as dark blue in Fig|3l 
The differences include also the variation in the position 
of Ed, thus, for completeness, we have recorded series of 
dl/dV (E,line) maps across the sample (see Fig. 0(a)- 
(d)). Both the intensity and the position of the observed 
local minimum depends on the spatial position on the 
investigated surface, which is very prominent on these 
plots. The signatures that seem to indicate Ed positions 
are placed between 0.25-0.55 eV (compare for example 
Figs 0(a) and (c)), in agreement with preliminary data, 
but also local variations in electronic structure can be 
easily noticed (slight deviations from definite values ob- 
served along the lines, for example the one in Fig|5](c)). 
Thus, the existence of large domains may be associated 
with different properties of visible crystallites, but the 
influence of other structural factors must also be taken 
into account. 

It seems that the doping level is determined by the 
combinations of many factors which can be identified de- 
pending on the scale of the observation. For example, the 
presence of repeated regions with well defined Ed above 
the Fermi-level and those revealing featureless spectra 
seem to be related mainly to structural properties of dif- 
ferent crystallites of the substrate (e.g. amorphous and 
crystalline domains), but also to its corrugation as well 
as to the local variations in the relaxation of the sur- 
face region. It seems, however, that in the case of the 
local changes illustrated in Fig. [5] rather the influence of 
last two factors should be dominant. In order to demon- 
strate this hypothesis, we have studied pure Au substrate 
to provide information concerning its structural and elec- 



tronic heterogeneity needed in further DFT simulations. 



C. STM/STS of pure Au substrate 

The STM/STS studies of gold have shown strong het- 
erogeneity of electronic properties which seems to be re- 
lated to the structural diversity as presented in Fig. [6^a). 
We have distinguished two main types of regions different 
in terms of LDOS - such as those illustrated in Fig. [6](b) 
and (c). Figure [6] (b) presents rather featureless spectra 
and reveals asymmetry between occupied and unoccu- 
pied states i.e. the LDOS is much higher at occupied 
states for all energies. This behavior is rather typical for 
amorphous gold substrates. In contrast, Fig. [6](c) shows 
spectra with continuous increase of LDOS for occupied 
states accompanied with clearly visible maxima located 
close to 0.4 eV below the Fermi level. This behavior seems 
to be the fingerprint of the Schockley surface state (SS) 
lying in the projected bulk band gap in T-L direction and 
observed on the Au(lll) surfac o 33 i 34 . It leads to the con- 
clusion that the gold substrate is partially covered by the 
crystallites with (111) orientation. 



IV. THEORETICAL MODELS OF 
GRAPHENE/GOLD INTERFACE 

The experimental studies of the sample revealed het- 
erogeneity in the electronic properties of graphene which 
seems to be domains of p-type and zero doping as well 
as the slight spread in the values of p-type doping within 
the respective regions. The STS spectra recorded in cer- 
tain regions of uncovered gold showed surface states typi- 
cal for (111) orientation, while featureless characteristics 
indicating the amorphous metal have been detected in 
other parts of the sample. Thus, the presence of do- 
mains can be associated with the structural diversity of 
the substrate. 

The interface between the graphene and a fee metal is 
usually modeled theoretically as a sheet of graphene ly- 
ing on a (111) metal surface. The (111) face has been 
elaborated due to the stability, symmetry and impor- 
tance of this configuration for practical applications^ —. 
First, we have performed the typical DFT simulations for 
graphene interacting with Au(lll). The considered sys- 
tem is presented in Fig. 1, while all the calculations' 
details were described in Sec. IIA. The calculations with 
semiempirical van der Waals corrections have shown that 
the Dirac cones are preserved in the electronic structure, 
but they are shifted of 0.03 eV above the Fermi level 
(AE F = E D - E F = +0.03 eV). The obtained band- 
structure is shown in Fig. [7] (b). One can notice that 
it agrees well with previous theoretical result G 13 i 17 . It 
should be also reminded that the level of doping is very 
sensitive to the distance between the substrate and the 
graphene sheets. Even its small increase induces higher 
values of doping. 
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FIG. 7: (Color online) (a) The electronic structure of free- 
standing graphene in 2x2 unit cell. The inset shows a Bril- 
louin zone with the marked T, K, and M high-symmetry 
points, (b) The band-structure of graphene on Au(lll). (c) 
The simulated LDOS spectra of a graphene/Au(lll) system 
for different positions above the surface. Fermi level lies at 
the energy value equal to zero. 



Figure [71 (c) presents the corresponding LDOS spec- 
tra obtained according to the method described in Sec. 
II A. Two curves have been simulated for different po- 
sitions above the surface. Their overall shape is similar 




FIG. 8: (Color online) (a) Top view of adsorption geome- 
try of graphene/ Au system with lateral rearrangement of Au 
atoms in the topmost substrate's layer, (b) Side view of 
graphene/ Au configuration with vertical rearrangements of 
Au atoms in the topmost substrate's layer. Carbon atoms are 
denoted as red (darker) balls, gold atoms as yellow (lighter) 
ones. 



to those obtained for graphene interacting with Cu(lll) 
and Pt(lll) surfaces 18 , but it seems that in the present 
case the features originating from complicated electronic 
structure of gold are better pronounced (compare Figs|6] 
(c) and [3(c)). Even in the LDOS profile calculated just 
above the surface (dashed line in Fig. [71 (c)) the states 
of Au(lll) can be visible in spite of the fact that carbon 
states dominate. The position of Dirac point (+0.03 eV) 
can be associated with the wide global minimum in the 
profile (dashed line), which is also confirmed by the band- 
structure plotted for this system (Fig. [71 (b)). These 
spectra seem to reflect the properties of graphene lying 
in the regions denoted as blue patches in Figj3j whose 
dl/dV profiles are presented in Figs. HI (a), (b) as green 
dashed lines. The STS measurements have shown the 
global minimum localized at the Fermi level which pro- 
vides a zero value of doping, in satisfactory agreement 
with theoretical results (+0.03 eV). The simulated pro- 
file (solid purple line in Fig (71(b)) is quite similar to the 
recorded one, but it should be noted that the Au(lll) 
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FIG. 9: (Color online) (a), (b) The band-structures of 
graphene/Au systems with vertical rearrangements of Au 
atomic positions. (c),(d) The corresponding simulated LDOS 
profiles. The zero energy is at a Fermi level. 



surface state is hardly detected by the STS in the regions 
covered with graphene (see Fig. 2]). 

Now, we introduce the interfaces containing strongly 
rearranged surfaces which can represent amorphous re- 
gions in the sample. We define the amorphous-type 
metallic surface as a set of atoms arranged randomly at 
the reasonable volume and relaxed. Such choice of the 
structure is the first step in the method of an ab initio 
random structure searching of the most stable configura- 
tions in the given systems—. It has been demonstrated 
that the relaxations of randomly distributed atoms pro- 



vide the most stable structure as the one the most fre- 
quently obtained among many trials. The other results 
of optimizations represent possible less stable configu- 
rations which should contain also amorphous-like sys- 
tems. Here, we consider only slight deviations from the 
(111) arrangement of atoms in a topmost layer of metal. 
This approach seems to be reasonable due to the ex- 
treme sensitivity of the doping level of graphene to any 
changes in the positions of atoms. The effect of such 
small variation should be then easily noticed. Moreover, 
we can distinguish the shifts of atoms with respect to the 
positions in (111) surface into lateral and vertical direc- 
tions. 

The geometry of the significant lateral rearrangement 
of atoms in a first layer is shown in Figj8] (a), whereas 
the slight vertical perturbation is illustrated in Fig|Sl(b). 
DFT calculations performed for the first system have 
shown that the lateral rearrangements do not lead to any 
variation in adsorption distance and to only a very slight 
change of the doping level (Ai?p=+0.08 eV) which might 
be caused by the local increase in distances between par- 
ticular C-Au pairs (C atoms are not longer on top of 
Au atoms). The band-structure corresponding to the 
configuration shown in Fig. [8] (b) is illustrated in Figj9] 
(a). One can observe that the conical point is shifted to 
AEp =0.28 eV, a value much larger comparing to results 
for Au(lll) surface (+0.03 eV). The Dirac point position 
can be also easily recognized in simulated LDOS profile 
(Figl9] (c)). In this case the shape of the spectrum is 
more similar to the experimental profiles which seem to 
correspond to the amorphous regions (solid red lines in 
Figsgl(a), (b)). 

The difference in doping between graphene interacting 
with Au(lll) and with amorphous-like surface presented 
in FigJH] (b) can be explained by the arguments of geom- 
etry and symmetry. The gold atoms in the first layer by 
definition does not lie in the same plane, because their 
positions are generated randomly in a given region. We 
can always distinguish the topmost atom, the middle and 
the lowest one. The distance between graphene and the 
plane of the topmost atom (d t0 p) is typically smaller (2.5- 
3.1 A) than the one optimized on perfect Au(lll) sub- 
strate (3.23 A), but the distance between graphene and 
the lowest atom is defined to be smaller than the aver- 
age one <i avg (3. 2-3. 7 A). Such configuration always leads 
to increase in doping, which could reach over +0.3 eV, 
but the accurate value strongly depends on the strength 
(amplitude) of the assumed perturbation. 

It should also be stressed that the average distance 
between the two topmost layers of gold atoms (d re i) 
usually changes comparing with the bulk layer spacing 
(^int=2.36A) and with the Au(lll) configuration (ap- 
proximately 2.5 A, see Table II in Ref.17). This effect 
is equivalent to average contraction or expansion of the 
layer spacing (the latter case is illustrated in Fig. [51(b)), 
which was demonstrated to strongly influence the doping 
levelii. In particular the increase in surface region inter- 
layer distances raises the doping, while the contraction 
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lowers its value (see Fig. 6 in Ref.17). It seems to be 
the explanation of the local modulations in the value of 
p-type doping observed in the spectra presented in FigJSJ 
but can enhance also the differences observed between 
the values in doping in separate regions on the sample 
(compare Figs. [5] (a) and (c)). 

It is worthwhile to note that the doping level of per- 
turbed system similar to the one presented in Fig|8] (b) 
does not achieve the maximal value of 0.55 eV observed in 
the experiment. It is due to the reorganization of bands 
in the vicinity of the K-point of the Brillouin zone which 
can be easily observed in Fig. [9] (a). Values of doping 
higher than +0.3 eV can be achieved for only slightly per- 
turbed systems, but having expanded two topmost layers 
of the gold. Such system behaves similarily to the one 
with increased interlayer spacing which was described in 
details in Ref.17 and illustrated in Fig. 6 therein. In Fig. 
[9] (b) and (d) we show, respectively, the calculated band- 
structure and LDOS profile for such configuration. The 
doping level achieves 0.41 eV, which is closer to the max- 
imal measured value. 

Finally, we would like to stress the role of the detected 
corrugation and the fact that the graphene does not map 
exactly the structure of gold. The configuration with 
graphene suspended on Au(lll) crystallites and above 
the amorphous regions provides the differences in the 
graphene-substrate distance in a natural way. Since the 
doping level is very sensitive to the graphene-substrate 
distance^, all the described variations in the Fermi level 
shift might be enhanced. 

V. FINAL REMARKS AND PERSPECTIVES 

To successfully employ graphene in the electronic de- 
vices and to construct systems promising for applications, 
it is important to understand the graphene-metal inter- 
action. The preliminary theoretical results 1 ^ showed the 
sensitivity of the doping of graphene to the geometry of 
the metallic substrate, which reveals potentially high tun- 
ability of such systems. Thus, we have used scanning tun- 
neling spectroscopy to study the properties of graphene 
deposited on conductive heterogenous substrate and per- 
formed DFT calculations to explain the obtained results. 

We have observed the domains of zero and p-type dop- 
ing across the graphene/ Au sample, but it should be 
stressed that the value of p-type doping differs between 



individual regions as well as varies slightly within each 
specific domain. Additional STS studies of pure Au sub- 
strate suggested that the existence of large domains is re- 
lated to the presence of both amorphous and crystalline 
regions with (111) orientation. Our theoretical studies 
seem to confirm this hypothesis. Moreover, even the 
changes in Au atomic positions within the topmost layer 
of the substrate lead to the transition from zero doping 
typical for perfect Au(lll) substrate to the Fermi level 
shift below the graphene's Dirac points. The distinct con- 
figurations of atoms in amorphous-type system provide 
different values of the p-type doping, in accordance with 
experimental data (compare Figs. [5] (a) and (c)). This 
factor seems to be also responsible for slight variation in 
the doping level within selected domains in the sample. 
The observed modulations might be associated with local 
changes in the average interlayer spacing in metal. Obvi- 
ously, the domination of every particular factor depends 
on the scale of the observation, thus the influence of those 
with the smallest impact can be detected only by studies 
of the internal electronic structure of the specific domain 
on the sample. 

Finally, we have shown that gold/graphene interface 
has no negligible influence on the electronic properties of 
graphene layer. This leads to the conclusion that this 
influence must be taken into consideration when the real 
graphene-based nanoelectronics devices are designed and 
built. 
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